Astronomy & Astrophysics manuscript no. half 


© ESO 2009 


December 16, 2009 





Letter to the Editor 



o 
o 

(N 
o 

Q 



The H-alpha luminosity function at redshift 2.2 * 

A new determination using VLT/HAWK-I 

Matthew Hayes 1 , Daniel Schaerer 1 ' 2 , and Goran Ostlin 3 

1 Geneva Observatory, University of Geneva, 51 chemin des Maillettes, 1290 Sauverny, Switzerland 
e-mail: matthew . hayes@unige . ch 

2 Laboratoire d'Astrophysique de Toulouse-Tarbes, Universite de Toulouse, CNRS, 14 Avenue E. Belin, 31400 Toulouse, France 

3 Oskar Klein Center for Cosmoparticle physics, Department of Astronomy, Stockholm University, 10691 Stockholm, Sweden 



Received September 01, 2009; accepted December 07, 2009 



ABSTRACT 



o 

u 

43 

9 L1, 
6 



> 

(N 

m 

On 
O 



We aim to place new, strengthened constraints on the luminosity function (LF) of H-alpha (Ho-) emitting galaxies at redshift z ~ 2.2, 
and to further constrain the instantaneous star-formation rate density of the universe (p*). We have used the new HAWK-I instrument 
at ESO-VLT to obtain extremely deep narrow-band (line; NB2090) and broad-band (continuum; A" s ) imaging observations. The target 
field is in the GOODS-South, providing us with a rich multi-wavelength auxiliary data set, which we utilise for redshift confirmation 
and to estimate dust content. We use this new data to measure the faint-end slope (a) of LFQia) with unprecedented precision. The 
data are well fit by a Schechter function and also a single power-law, yielding a = (-1.72±0.20) and (-1.77±0.21), respectively. Thus 
we are able to confirm the steepening of a from low- to high-z predicted by a number of authors and observed at other wavelengths. 
We combine our LF data-points with those from a much shallower but wider survey at z ~ 2.2 (Geach et al. 2008), constructing a 
LF spanning a factor of 50 in luminosity. Re-fitting the Schechter parameters, we obtain logL* = (43.07 ± 0.22) erg s , log0, = 
(-3.45 ± 0.52) Mpc~ 3 ; a = (-1.60 ± 0.15). We integrate over LF(Ha) and apply a correction for dust attenuation to determine the 
instantaneous cosmic star-formation rate density at z ~ 2 without assuming a or extrapolating it from lower-z. Our measurement of 
Pi, is (0.215 ± 0.090) Mq yr" 1 Mpc~ 3 , integrated over a range of 37 < log(L Htt /erg s" 1 ) < 47. 
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1. Introduction 

Since th e evolution in the rate of cosmic star -formation was first 
plotted dLiliv et alj|1996t iMadau et al"1ll996l) . the literature has 
been awash with studies adding further points to the diagram. 
To estimate this quantity, one typically needs to find a num- 
ber of galaxies by a given method and estimate their density 
in both space and luminosity (i.e. the luminosity function, LF). 
Integration of LF ■ L - dL thus provides the volume-averaged lu- 
minosity density (pi) and, if L is a suitably calibrated measure 
of star-formation rate (SFR), pi converts directly into the volu- 
metric rate of star-formation (p*). 

A favourite tracer among low-z observers for estimating star- 
formation is the Ha emission line due to its simple physics, high 
intrinsic brightness, and convenient rest-wavelength at 6563A. 
At z > 2.5, however, the Ha line shifts out of the /T-band, mak- 
ing it highly inefficient for galaxy evolution studies. Fortunately 
at this z, selection by either the Lya line or continuum (BM/BX, 
BzK) criteria becomes possible, but unfortunately a different 
population of galaxies may be recovered, introducing biases 
in selected star-formation rate/history, dust content, etc. More 
specifically with regard to SFR, continuum luminosities come 
to equilibrium (and are therefore calibrated) over very different 
time-scales to nebular lines: > 100 and ~ 10 Myr respectively, 
impl ying that lines a re more sensitive to the instantaneous SFR 
(e . g. iKennicuttl 1 998b . Lya enables surveys to go much further in 



redshift and has identical production physics, but is a resonance 
line and undergoes a complex radiation transport, which renders 
it an unreliable tracer of intrinsic properties. From a purely phys- 
ical perspective, Ha is a far preferable tool. 

Luminosity functions are typically parameterised and com- 
pared using the Schechter functionj in which the luminosity 
distribution below (above) the characteristic luminosity, L*, is 
dominated by a power-law (exponential). When fitting, the three 
Schechter parameters are degenerate, and strong constraints are 
only obtained by sampling above and significantly below L*. 
Many attempts have been made to pin d own LF(Ha) and it s 
evolution with redshift from z ~ (e.g . iGallego et al.l|1995[ 
through inte rmediate-z ( Yan et al. 19991; iHopkins et al.l |200< 



iTresse et dj2002t ISobral et al.l2009tlShim et al.l2009l) . with the 
first z — 2 limits placed by Bun ker et al.l d 19951) . It is only in 
recent years that z — 2 surveys have been fruitful, with the 
most significant study being that of lGeach et al.l d2008l hereafter 
G08). However, this study was a wide, shallow survey designed 
to find the bright objects and determine L+, but does not permit 
an estimate of the faint-end slope. 

With the wide-field (7:5x 7 :5) HAWK-I instrument 
dPirard et alj|2004l ICasali et al.ll2006l) at ESO-VLT we have ob- 
tained the deepest narrowband Ha observations to date as part 
of a programme to study Ha and Lya emitting galaxies from 
a single volume at z ~ 2.2 (Hayes et al., submitted). These Ha 
data alone enable us to study the faint-end of LF(Ha) at unprece- 



* Based on observations made with ESO Telescopes at the Paranal 
Observatory under programme ID 081.A-0932 
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dented depths, tightening constraints on the overall LF (Ha), and 
providing the content of this Letter. In Sect. [2] we briefly de- 
scribe the data, reduction, and selection; in Sect. [3] we present 
the z ~ 2 Ha luminosity functions we derive and the constraints 
on pi,; and in Sect.|4]we briefly summarise. Throughout we adopt 
a cosmology of (Z/q,Q m ,Qa ) = (70 km s" 1 Mpc" 1 , 0.3,0.7), 



the SFR(Ha) calibration of iKennicuttl d!998l) . derived for a 
Salpeter initial mass function in the mass range (M\ a ,Mhi) = 
(0.1, 100)Mq and solar a bundances, and the AB magnitude sys- 
tem dOke&Gunnlll983l) . 

2. The data 

2.1. Observations and reduction 



A field in the GOODS South dGiavahsco etalj|2004 was se- 
lected to maximise the depth and quality of the auxiliary data. 
The central pointing was a =03:32:32.88; 5 =-27:47:16; p.a.= 
-44 deg. The field was observed in service mode between 08 
September and 21 August 2008 on a single dithered pointing. 
The NB2090 filter (A c = 2.095/im; AA = 0.019/mi) was used 
to isolate emission line candidates, using a total integration time 
of 60000 seconds. For the continuum we obtained 7 500 s in 
HAWK-I K s band, which we combined with the publicly avail- 
able ISAAC K s data from the GOODS/ESO Imaging Survey. 

Full details of the reduction process are beyond the scope of 
this Letter and will be presented in a forthcoming article (Hayes 
et al. 2010). For brevity, the EsoRex HAWK-I pipeline was used 
on the individual frames for bias-subtraction, flat-fielding, and 
subtraction of the sky with temporally adjacent image pairs. 
Custom scripts were then used to mask cross-talk artifacts, reg- 
ister and co-add the individual frames. We modelled the point 
spread function (PSF) in the final image and determined a see- 
ing of 0.89". To estimate the limiting magnitude we added arti- 
ficial point-sources (full width at half maximum set to the mea- 
sured seeing) to the images with the addstar t ask in noao/iraf 
and te sted their recovery using SExtractor dBertin & Arnoutsl 
1996; see Sect. 12.2b . It should be noted that the limiting magni- 
tude for extended objects will be somewhat shallower, but also 
that at z = 2.2, one seeing disc corresponds to a physical scale 
of 7.4 kpc. Using SExtractor mag_auto magnitudes we de- 
termined a 5<x limiting magnitude of 24.6 in NB2090. By com- 
puting the product of this flux density and A/Ijvt#2090> this corre- 
sponds to a line flux of 6.85xl0 -18 erg s _1 , if all the NB2090 flux 
comes from Ha and falls at the peak of the filter throughput. This 
assumption, valid for perfect top-hat bandpasses, holds well for 
our bandpass, which has steep edges by narrowband standards: 
the full width at 80% transmission is over 80% the FWHM. This 
limiting flux corresponds to SFR=1.9M Q yr _1 . The K s 5<x limit 
is 25.4. 

2.2. Selection of H-alpha emitting candidate galaxies 

Source detection was performed in the on-line image using 
SExtractor, where we required a minimum of 12 contiguous 
pixels (plate-scale^/ 106 px _1 ) to reside above the background 
noise by a signal-to-noise {S /N = 3). Aperture-matched K s 
photometry was performed using "double image" mode. We se- 
lected emission line candidates based upon two criteria, the first 
of which was a minimum rest-frame equivalent width, EWn a - 
20A. Furthermore, in order to prevent over-contamination from 
,K" s -faint objects scattering into the colour-selection region, we 
required the narrowband flux to be a factor of 2 = 5 greater than 
the noise in the continuum image. Figure Q] shows the colour- 
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Fig. 1. The selection function for candidate Ha emitters. Black 
dots represent all the objects detected in the NB2090 image. The 
red dashed line shows the colour cut of EWu a = 20 A, and the 
green solid line the limit of Z = 5. The blue points show the 
selected candidate galaxies. The inset histogram shows the nar- 
rowband magnitude distribution of all candidate emission line 
galaxies in blue, and confirmed z s 2.2 Ha emitters after the 
rejection of interlopers (see text) in magenta. 



magnitude selection diagram, including the cuts in EW\\ a and 

After inspecting each candidate by eye, we found 152 objects 
that could potentially be emitters of any emission line between 
Paa at z ~ 0.1 1 and Lya at z ~ 16. We then cross-correlated our 
candidates with the merged z- an d A^-selected GOODS-MUSIC 
catalogue of ISantini et al.l d2009l) . finding 143 matches. The ob- 
jects for which counterparts were not detected all appear to be 
genuine narrowband-excess objects, but with equivalent width 
lower limits so high that their stellar continua are too faint to 
be detected in the ISAAC K s or ACS z-band images. GOODS- 
MUSIC provided us with spectral energy distributions (SEDs) 
between the [/-band and Sptizer-MlPS 24ytzm. To assemble the 
final catalogue we first examined the spectroscopic redshift mea- 
surements in the GOODS-MUSIC catalogue. If the spec-z had a 
quality flag of "very good" or "good" and was consistent with the 
redshift interval defined by the NB2090 bandpass (2.178 < z < 
2.207) the galaxy was included (1 galaxy). Using the same qual- 
ity criterion, we excluded objects if they fell outside that z range 
(26 galaxies). This ratio of 26 galaxies excluded versus one in- 
cluded may superficially appear to show the selection method in 
a rather negative light, but the bias can easily be explained by 
considering the selection methods of the many studies that fol- 
lowed up the GOODS imaging spectroscopically. We searched 
the spectroscopic catalogues for all galaxies in GOODS-MUSIC 
with spec-z in the range covered by FWHMnbxm) and spec-z 
flags of good or better, finding only a single object that could be 
detected by our survey: the one we do find, followed up from 
BzK selection. In contrast the catalogue is rife with spectra tar- 
geting the z = 0-2 domain and Lyman-break galaxies at z > 3.5, 
and the combined catalogue is heavily biased against z = 2 — 3. 
We find typically as interloping lines: z ~ 3.2 [Om]; low-z Paa, 
PaB and lines in the higher order Brackett series; and z ~ 0.28 
[Fe n] emitters, although a handful of more unknown or inter- 
esting lines are also recovered. The number of emission line 
candidates is of sufficient interest to warrant further examina- 
tion, which will be the topic of a forthcoming paper. For ob- 
jects with no spec-z or an uncertain flag, we used the Hyper- 
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Fig. 2. Left: Ha luminosity functions. The black points show the bins created from HAWK-I data obtained for our programme 
with error-bars derived from Poisson statistics and Monte Carlo incompleteness simulations. The best-fit Schechter function is 
shown by the black solid line, while the black dotted line represents th e best fitting simple power-law. Othe r lines show variou s 
lumi nosity function s from the lite rature, scale d to our cosm ology. H00=dHopkins et al]l200"o} all po i nts), G 08= lGeach et al.l ((2008), 
G95= lGallego etalj d 1995b . T02= lTresse etail d2002l) . So09= ISobral et alj d2009l) . Sh09= IShim eTaT] d2009t) . Dashed vs. solid repre- 
sentation shows the 3<x limiting survey depth. Right: The complete lum inosity funct i on co mbining these HAWK-I data with the 
points from G08. The shaded cyan area shows the LF(Ha) derived by iReddy et alJ d2008l) from the modelling of < z >= 2.3 
colour-selected galaxies. 



z photometric -redshift code (Bol zonellaetal.ll2000h. modified 
to inc lude nebular continuum and lines dSchaerer & de Barrosl 
2009). Here we selected objects that had lcr errors on their phot- 
z consistent with z = 2.19. Combined with the spectroscopically 
confirmed objects this gave us 55 z ~ 2.2 Ha-emitters. For se- 
curity, we tested our gal axies against the B zK colour criterion of 
(B-z)-(z-K)> -0.2 dDaddi et alj2004l) to insure that our ob- 
jects would be classified as star-forming objects at z > 1.4. Only 
four of our 55 objects do not satisfy this criterion, but when ex- 
amining the BzK colours of the entire narrowband-excess sam- 
ple (143) we found an additional 7 objects that do. We note also 
that the BzK criteria neglect the fact that the B and K photometry 
may be contaminated by the Lya and Ho- emission lines, respec- 
tively. Since the number of objects that are classified/declassified 
by BzK colours are (/.) so similar, (ii.) only a small (~ 10%) frac- 
tion of our sample, and (Hi.) possibly contaminated by lines, we 
opt to leave our phot-z selected sample unchanged. Finally, we 
cross-corre lated our sample against the 1 Mega-seco nd Chandra 
Catalogue dGiacconi et al.ll2002t iRosati et alj|2002l) . but no ob- 
jects in which the Ho- production is obviously dominated by 
an active nucleus were found, including the objects missed by 
broadband cross-correlation. Our final catalogue comprises 55 
objects. 

3. Results and Discussion 

3.1. The redshift 2.2 H-alpha luminosity function 

All Ha photometry is corrected for the contribution of 
[Nn]/U6548,6583A lines. In local starbursts the [Nn]/Ha ra- 
tio varies strongly with metallicity (Z) between ~ and 0.6, 
and [N n]/Ho' may plausibly be estimated from Lu a through the 
luminosity-metallicity relationship. However a significant offset 



Table 1. Luminosity function parameters 



Param. 


Unit 


HAWK-I all" HAWK-I a " 


Combined c 




erg s~ J 


43.22 + 0.56 


43.07 +0.22 




Mpc- 3 


-3.96 ±0.94 


-3.45 ± 0.52 


a 




-1.72 + 0.20 -1.77 + 0.21 


-1.60 + 0.15 



" 3 Schechter parameters: HAWK-I data only. 
b a only: HAWK-I data only. 

c 3 Schechter parameters: HAWK-I data + G08 points. 



in the L - Z relationship is found at z ~ 2 dErb et alj|20 06a). 
and an application of the l ocal relationship at z ~ 2 would be 
insecure. Furthermore, the lErb et all d2006al) study is based on 
galaxies an order of magnitude brighter than those found here, 
and we dare not extrapolate to these luminosities. Thus we ap- 
ply a conservative single correction of [N n]/Ha= 0. 1 1 7, derived 
by ave raging over the faintest four bins in the data of Er b et al.l 
d2006ah . 

LF(Ha) is created by binning all selected objects by lumi- 
nosity, with errors derived from the standard Poisson statistic and 
Monte Carlo simulations of incompleteness. We fit the Schechter 
function using a standard^ 2 minimiser, determining errors from 
1 000 re-drawn Monte Carlo realisations. The best-fit Schechter 
parameters are presented in Table Q] and are shown graphically 
with other determinations from the literature in Fig. [2] 

The brightest galaxy in our sample has a luminosity of 
3.5 x 10 42 erg s~' , which is a factor of 4.5 fainter than the value of 
Li, we derive, and even allowing for the error on L*, we clearly 
are not probing this luminosity domain. This galaxy is a factor 
of 2 fainter than the L* = 6.7 x 10 42 erg s^ 1 derived by G08 
and demonstrates the inadequacy of such a small survey to con- 
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Fig. 3. The cosmic star-formation rate density. Open grey pen- 
tagons show a ll of the dust-corrected points compiled by 
Hopkins (2004) that were derived from anything oth er than Ha. 
Filled red squares show dust-corrected points from Shi m et all 
(120091) . all of which are Ha-based. The black point is the mea- 
surement from this study and the open blue point is the result 
from G08, artificially shifted to higher-z to allow it to be distin- 
guished from our own. 



strain all Schechter parameters. Indeed, the statistical error on 
the degenerate parameter 0+ spans an entire dex and it is appar- 
ent that the only parameter we can reliably constrain is a, since 
all our Ho- emitters have luminosities that place them firmly in 
the power-law region of the LF . For security, we proceed to ex- 
amine the faint-end only by fitting of a simple power-law, using 
exactly the same Monte Carlo realisation data as previously. Our 
derived faint-end slope is found to be -1.77+0.21, also presented 
in Table [1] and exhibits a near-identical mean and errorbar. 

Some authors find a steepening in the faint end slope of 
LF(Ha) moving from z = to z > 1 with similar results fou nd in 
UV-sel ected catalogues at z = 2 to 3 dReddy & Steidell2009l) and 
z = 6 dBouwens et alj|2007l) . In addition. iHopkins et al.l T2000) 
find a steepening in a for Ha emitters at z ~ 0.7 to 1 .8 by slitless 
spectroscopy with the HST/NICMOS, with a = -1.60 or 1.86 
depen ding on the selection criteria. However, iTrenti & Stiavellil 
(2008) note that for such small field-of-view studies, an arti- 
ficial steepening of the faint-end slope may manifest itself if 
under-dense regions are targeted. Examining the redshift distri- 
bution of galaxies in the whole GOODS field using both spec-z 
and phot-z reveals z = 2.2 to correspond to neither a signifi- 
cant over- nor under-density. After deriving a - -1.77 + 0.21, 
we are able to confirm this increase in the relative abundance 
of lower-luminosity galaxies at epochs of peak cosmic star- 
fo rmation. This confirm ation comes at the 3<x level with respect 
to lGallego etaild 1995b (z = 0), who find a = -1.35, and our 
results are c onsistent with no e volution within lcr to z = 0.84 
(a = -1.65:ISobral etalJ l2009h and 0.8 < z < 1.7 (a = -1.86; 
HopkinseLalEoOO). 



3.2. Combined LF(Ha) and the star-formation rate density 

As stated above, our data are not sufficient to constrain the 
brighter end of LF(Ha). To address this we take the LF (Ha) 
data-points from G08 and combine them with our own. We ex- 
tract the LF data-points, which have been corrected for the [N n] 



using the z ~ estimate of Kennicutt & Kent ( 1983), and reap- 
ply the z ~ 2 correction of Erb et al.l d2006al) for consistency with 
our own points. The resulting LF can be seen in the right panel 
of Fig. [2] with the best-fitting Schechter parameters, derived fol- 
lowing exactly the same method as in Sect. 13.11 in Table Q] We 
also show an estimate of th e observed (i.e. du st-uncorrected) 
LF(Ua) at < z >= 2.3 from lReddv etalj d2008l) . inferred from 
colour-selected samples by equating dust-corrected SFR(UV) 
with SFR(Ha) and re-applying the effects of dust attenuation to 
Ha. The observed and inferred LF(Ha) agree well in general, 
particularly around L*, with only a mild deviation at the faint 
and (unsampled) bright ends. 

In order to estimate the cosmic star-formation rate density, 
we need to correct our Ha fluxes for dust attenuation. To do this 
we re-fit the SEDs of all the Ha-emitters, fixing their redshift to 
2.19. Ideally we would correct our luminosity function bin by 
bin, but no obvious correlation is observed between Lu a and Ay. 
Thus we adopt the average Ay f or the sample of 1.19 + 0.55, 
assuming the attenuation law of ICalzetti et al.l (120001) . We do 
not employ a correction factor for the effects of diffe rential at- 
tenuatio n of stellar and nebular light, since at z ~ 2 lErb et al] 
(l2006bl) find SFR(UV) equivalent to SFR(Ha) without this fac- 
tor. Ay- 1.19 equates to A 6563 ^ = 0.977 for the fainter bins in 
LF(Ha), identical to the assumption of 1 mag used by G08 at the 
brighter end, and we corr ect for 1 mag of a ttenuation through- 
out. It is noteworthy that iReddv et al .1 (120081) find a maximum 
likelihood E B - V of 0.12 (A 6563A = 0.4, Calzetti dust) for BX- 
selected galaxies at < z >= 2.3. The significantly different ex- 
tinction determined from Ha selection may be due to the very 
different selection functions: indeed EWho is independent of 
dust extinction, whereas BX selection cuts the multi-broadband 
sample based on U - G colours, and thus it is not unreason- 
able to infer that galaxies with very different extinctions sur- 
vive the respective cuts. The errors derived on the Schechter 
parameters are strongly correlated, and thus we compute p + 
from each of the raw Monte Carlo realisat ions previously de- 
scribed. For consistency with lHopkinsI d2004l) we integrate over 
the luminosity range 37 < log(LHo-/erg s -1 ) < 47, where we 
find a value of p L = (2.72 + 1.14) x 10 40 erg s _1 Mpc~ 3 and 
p* = (0.215 + 0.090)M© yr" 1 Mpc~ 3 . This however repre- 
sents a significant extrapolation over luminosity, and we also 
present the same integration performed over the range covered 
by our bins in total LF(Ha): 41.5 < log(L H »/erg s -1 ) < 43.1. 
This value, much less dependent on the Schechter parameter fit- 
ting, results in p L = (1.76 + 1.02) x 10 40 erg s _1 Mpc~ 3 and 
p* = (0.139 + 0.081) M Q yr" 1 Mpc~ 3 . We show the result 
in Fig. [3] along with two compilations of p* from the litera- 
ture. We in c lude t he dust-corrected Ha survey da t a co mpiled by 
Shi m et al .1 (120091) . and points from the iHopkinsI (|2004|) compi- 
lation based upon many different observations where we have 
excluded all determinations based upon Ha. We also indicate 
separately the point of G08, which has been artificially moved 
in redshift to separate it from our own. Our point agrees well 
with that of G08, which is initially surprising, given that we ob- 
tain a significantly steeper faint-end slope and marginally higher 
L*. The culprit behind the similar p+ yet differing (L+,a) is nat- 
urally the overall normalisation of the luminosity function, 0*, 
which we calculate to be a factor of 4 lower. Furthermore the er- 
rorbar produced by our analysis is actually larger than that pub- 
lished by G08, despite our LF extending an order of magnitude 
fainter in luminosity. These errorbars, however, are statistical er- 
rors that result from averaging over many Monte Carlo realisa- 
tions, and our fitting engine includes an additional parameter (a) 
- fixing one parameter introduces a systematic error that is not 
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accounted for by Monte Carlo. Had we locked a to its best-fit 
value and fitted L* and 0+as usual, our statistical errorbar on 
p* would be a factor of 3 smaller. Finally it is possible that our 
selection criterion of EWu a > 20A could cause us to underes- 
timate the total luminosity density. To assess the impact of this 
we examine the EWu a distribution of our sample and fit an expo- 
nential of the form oc exp(EW/EW Q ), finding EW Hafl = 1 15A. 
Assuming the two extremes of EWn a to be independent of Lu a 
and also of the continuum luminosity, we find that our selection 
criterion misses between 1.3 and 16% of the integrated luminos- 
ity. We also examine the EWh<i distribution of nearby galaxies 
in the Sloan Digital Sky Survey, selecting the complete spectro- 
scopic sample at z < 0.1, and find that removal of galaxies with 
EW Ha < 20A cuts 19% of the Ha luminosity density. 



4. Summary 

We have used the new HAWK- 1 instrument mounted at ESO-VLT 
to carry out a blind, narrowband survey for Ha emitting galaxies 
at a redshift of 2.2. This is the deepest unbiased survey carried 
out at this redshift to date and enables us to estimate the faint- 
end of the Ha luminosity function, a parameter hitherto merely 
assumed at this cosmic epoch. The target is the GOODS-South 
field, offering us a rich, deep multi-wavelength ancillary data- 
set, in which we find 55 Ha emitters. 

1. We construct a luminosity function from the HAWK-I data 
which we find is well fit by a Schechter function with pa- 
rameters of logL* = (43.22 + 0.56) erg s"\ log^* = 
(-3.96 + 0.94) Mpc~ 3 , a = (-1.72 ± 0.20). Fitting a sin- 
gle power-law component to all the LF points gives a very 
similar value of a = (-1.77 + 0.21). This confirms the steep- 
ening of the faint-end of the galaxy luminosity function out 
to z = 2, which is predicted and observed at other wave- 
lengths. 

2. We combine our luminosity function bins with those from 
the much wider but shallower survey of G08, which yields 
the best-sampled Ha luminosity function at this redshift. 
This is also well described by a Schechter function with the 
parameters of logL* = (43.07 + 0.22) erg s -1 , log0* = 
(-3.45 ± 0.52) Mpc~ 3 , a = (-1.60 ± 0.15). 

3. We apply corrections for dust attenuation derived from mod- 
elling of the stellar spectral energy distributions and in- 
tegrate the combined luminosity function to obtain an in- 
stantaneous cosmic star-formation rate density of (0.215 + 
0.090) Mq yr _1 Mpc~ 3 . This provides the most robust 
emission-line estimate to date at this redshift. 
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